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A priori testing of subgrid-scale 
models for the velocity-pressure 
and vorticity-velocity formulations 

By G. S. Winckelmans 1 , T. S. Lund 2 , D. Carati 3 AND A. A. Wray 4 


Subgrid-scale models for large eddy simulation (LES) in both the velocity-pressure 
and the vorticity-velocity formulations were evaluated and compared in a priori 
tests using spectral Direct Numerical Simulation (DNS) databases of isotropic tur- 
bulence: 128 3 DNS of forced turbulence ( Re\ = 95.8) filtered, using the sharp cutoff 
filter, to both 32 3 and 16 3 synthetic LES fields; 512 3 DNS of decaying turbulence 
( Re\ = 63.5) filtered to both 64 3 and 32 3 LES fields. Gaussian and top-hat fil- 
ters were also used with the 128 3 database. Different LES models were evaluated 
for each formulation: eddy- viscosity models, hyper eddy- viscosity models, mixed 
models, and scale-similarity models. Correlations between exact versus modeled 
subgrid-scale quantities were measured at three levels: tensor (traceless), vector 
(solenoidal 'force’), and scalar (dissipation) levels, and for both cases of uniform 
and variable coefficient(s). Different choices for the 1/T scaling appearing in the 
eddy-viscosity were also evaluated. It was found that the models for the vorticity- 
velocity formulation produce higher correlations with the filtered DNS data than 
their counterpart in the velocity-pressure formulation. It was also found that the 
hyper eddy- viscosity model performs better than the eddy viscosity model, in both 
formulations. 


1. Velocity-pressure formulation and models investigated 

Consider the Navier- Stokes equations for incompressible fluid in the velocity- 
pressure formulation: 


du x d dP d du t 

~dt dxj ( U,Uj) + dTi~ "dxjdxj ' 

Filtering, using a low-pass filter G of characteristic length A, 

V-’(*) = J 4>(v)G ip'( x ) = -4’{x) , 
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with %j> the filtered value and if?' the remainder, leads to the following evolution 
equation for the filtered velocity field: 

du t d , x dP d dUi 

w + er, + = 

which is rewritten either as 


dui d dP dr tJ d du, 

~m + dT, + + = ' 


or as 


dui 


dP df 


'] 


d dui 


~5t + SJ (u, “ ;) + S7 + &7 = "dx, Sx~ 

where the ‘subgrid-scale stress 1 (sgs) tensor is defined as 


UiUj — C jj -f- Riy 4* i 
Uji/j ^ C^ij 4“ Rjj ) 

with the usual definitions for the cross term, the Reynolds term, and the Leonard’s 

term: 

Cij = Uiu'- + u\uj , 

= u'u' , 

L-j = UjUj = — (u t i7j ) l . 

Here, Cjj and R l} are purposely written with the ‘overline’ as they are the filtered 
value of some quantity. L\> is purposely written with the ‘prime’ as it is the remain- 
der, after filtering, of some quantity and thus contains high spatial frequencies. 

The notation r,j is somewhat misleading (but will nevertheless be retained). 
Indeed, although it is the filtered value of some quantity, it is not necessarily the 
result of filtering r tJ . In the case of sharp cutoff filter in wave space, 4 y = V’ and 
t/>' = 0 so that Tij is indeed the result of filtering r x j. In the case of smooth filters 

such as the Gaussian, the filtering of r %) produces C tJ + R tJ + 17 ij, which is not 
equal to C tJ + R tJ . 

When doing LES with a computational grid which is of the same size as the 
assumed filter size A, one cannot accurately evaluate quadratic terms such as (uiilj ) 
using the coarse grid only, as such quantities have a high frequency content. One 
can only hope to resolve, on the coarse grid, quantities such as (and even that 

requires using an appropriate reconstruction scheme, e.g., the need for dealiasing 
in spectral codes). The second form of the filtered equation is thus the one to 
consider in LES computations. It is assumed throughout this paper unless otherwise 
specified. (Notice that, if one were to use a computational grid smaller than the 
assumed filter size, then quantities such as TtiUj could be partially resolved.) 



def 


T -J 


UiUj — 


def 


¥ ij 


UiUj — 



A priori testing of sgs models in the u — p and u; — a formulations 


311 


The trace of the subgrid-scale tensor does not influence the dynamics of the 
filtered flow and is usually lumped into the pressure term: 

T* } = Tij - i T kk Sij , P* = P + ^ T kk , 

du { d __ dP* dr*, d dui 

-dT + dTj + + ^7 = v d^dTi ' 

Only the subgrid scale force, the divergence of the subgrid-scale tensor /* = 
needs to be modeled: 3 degrees of freedom instead of 5, (or 6 if one were also 

interested in modeling r**). The modified pressure is solution of 

a &T d d df* 

dx t dx t dxi dxj U,Ujl + dx{ 

Finally, the solenoidal (i.e., divergence-free) part of the subgrid-scale force is the 
only one that affects the flow dynamics. Defining g* as the solenoidal part of /*, the 
other part being the gradient of a potential <f> which is lumped into a new ‘pressure 1 , 
V — P + <f>, we write: 


du x d dV* d du l 

~dT + d^ {u,Ui) + ~d^ +9 ' -"^ 7 ^ 7 ’ 
d &P* _ d d 
dx x dx x dxi dxj 

Correlations of different LES models with filtered DNS data in isotropic turbu- 
lence were obtained and investigated. This was done at three levels: tensor level 
(traceless sgs tensor), vector level (solenoidal force), and scalar level (dissipation), 
for two different DNS data sets, using the sharp cutoff filter in wave space (with 
spherical truncation): 

a) 128 3 DNS of forced isotropic turbulence (Ke\ =. 95.8) that was filtered to both 
32 3 and 16 3 synthetic LES fields, see Fig. 1 . 

b) 512 3 DNS of decaying isotropic turbulence (Rex = 63.5) that was filtered to both 
64 3 and 32 3 synthetic LES fields, see Fig. 2 . 

In addition, correlations at the tensor level were also obtained when using smooth 
filters (here applied in wave space) with G(x/ A) = n*=i 3 G(k A) = 

FUi 3 G(ki A), such as the Gaussian (of same standard deviation as the top-hat), 



and the top-hat, 



G (ki A) = exp [ — 


(-HW)' 


*x_ 

A 



0 otherwise ; 


G(k,A) 
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Figure 1 . Spectrum of the 128 3 DNS {Rex = 95.8). Also shown is the cut used 
to produce the 32 3 synthetic LES field. , E(fc); , A'“ 5 ^ 3 . 


to generate synthetic LES fields from the 128 3 DNS data (i.e., DNS with t max = 64). 
The cutoff wavenumber was set to A* ma x = 16, hence A = 7r/fc m ax = tt/ 16. Notice 
that, with smooth filters, the synthetic LES field still contains contributions from 
all original modes. 

Many different LES models were investigated: 

Model 1 (eddy-viscosity type, tensor modeling): 


T*ij M — —2 Ft Sij 


- 1 (dxii du,\ 

Wlth S ^ = 2(^7 + ^J 


Different choices for the 1/T scaling that appears in the eddy- viscosity were 
investigated: 


v% 

- C A 2 


(a) 

Vt 

= C A 2 


(*) 

V% 

< 

O 

II 

('/A*)' 1 ’ , 

(c) 

Vt 

= C A 2 

(2u>iSijiJj) /(uJ.u?,) , 

(d) 

Vt 

< 

O 

II 

(2 uJjSijuJj) I/3 , 

(e) 


where e is the rate of energy transfer within the inertial range (assumed constant) 
and u) = V x u is the large-scale vorticity. The first choice is the classical Smagorin- 
sky’s scaling based on local dissipation by the large scales (e.g., see reviews by 




A priori testing of sgs models in the u — p and u> — u formulations 


313 



Figure 2. Spectrum of the 512 3 DNS (Re* = 63.5). Also shown are the cuts 
used to produce the 64 3 and 32 3 synthetic LES fields. • : Comte-Bellot and Corrsin 
experimental data at Re* = 65.1. k normalized by Kolmogorov scale g. 


Rogallo & Moin 1984, Lesieur et al. 1995). Consider the eigenvalues Ai, A 2 and A 3 

of the rate-of-strain tensor, with A] + A 2 + A 3 = Skk = 0. The scaling then produces 

1 /2 

an eddy- viscosity proportional to (2 (Af + \\ + A3)) 

The second choice is based on local enstrophy of the large scales (e.g., Mansour 
et al 1978). Recalling the identity 

d d 

2 SijSij = “f* 2 ^ (Ujifj) , 


together with the Poisson equation for pressure, it appears that, to first order, the 
two scalings differ by the local pressure Laplacian (which can have either sign). 

The third choice was proposed by Carati et al (1995b) and is referred to as 
the ^Kolmogorov’ scaling. If, following Smagorinsky, local equilibrium between the 
rate of energy transfer within the inertial range is identified with the subgrid-scale 
dissipation e % — 5,y, one recovers the classical Smagorinsky ’s model for the 

eddy-viscosity. The Kolmogorov scaling has the practical advantage over the other 
models that fewer filtering operations are required when implementing the dynamic 
procedure as in Germano et al (1991), Ghosal et al (1992), Moin & Jimenez 
(1993), Moin et al (1994), Ghosal et al (1995), Carati et al (1995a). When 
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developing the dynamic procedure in the present formulation, one obtains 


-2A 4/3 


1/3 


C'S t1 -C'S t 


= Uj - u t Uj ^ 


where ip is an additional Test filter 1 of the LES field, at scale A, where C' stands 
for the dimensional product C e 1//3 , and where it has been assumed that both filters 
lie in the inertial range so that e is indeed the same at both filter sizes. Since 
the dynamic procedure also assumes that the model coefficient C is invariant with 
filtering scale in that range, it turns out that C' is also invariant with filtering scale 
and is, in fact, what is determined by the dynamic procedure. 

The fourth scaling was proposed by Winckelmans (1995) and is based on the 
relative rate of change of the large scale enstrophy due to 3-D stretching, 

Ut Uju 

This scaling selects the eigenvalues used to compute the eddy-viscosity according 
to the relative orientation between u> and the principal axes (eigenvectors) of the 
rate-of-strain tensor. Indeed, writing the components of uj in the system of prin- 
cipal axes as (a; i,u; 2 ,u> 3 )i this scaling produces an eddy- viscosity proportional to 
2 (AjaJi + A2 u ?2 + A3UJ3) / (u;? +W2 +^3). Hence a vorticity-weighted average of 
the eigenvalues is used to produce the eddy- viscosity. This scaling produces a neg- 
ative eddy viscosity scaling in regions where enstrophy is decreasing. To ensure 
positivity, one can use either |d| or = max(d, 0). 

For completeness, a fifth scaling is considered, which is based on the rate of change 
of enstrophy, hence an eddy viscosity proportional to (2 (AiuJ, + A 2 u7 2 T A3LJ3)) 1 ''' 3 . 

Models other than the classical LES model (Model 1) were also investigated: 
Model 2 (eddy- viscosity type, solenoidal force modeling): 

r M = - v x (*W) . 

Model 3 (hyper eddy- viscosity type, tensor modeling): 

r XJ M = 2A 2 f ( v 2 s 0 . 

Model 1+3 (mixed eddy- viscosity and hyper eddy viscosity type, tensor modeling): 

r*+ = -2 (v n Sij - A 2 u t2 V 2 S U ) . 

Model 1+4 (mixed, simplest non eddy-viscosity type, tensor modeling): 

T ij M — ~ 2 (vt Sij + C 2 A 2 (5,/ S/j) ) . 

We are unaware of any published results describing Model 3 or Model 1+3. How- 
ever, Cerutti and Meneveau (1996, private communication) have also considered this 
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model. Model 1+4 is one of the many models investigated by Lund and Novikov 
(1992) where they considered all possible models with the sgs tensor function only 
of the strain and rotation rate tensors. 

Notice that one can even build models that give a decent (i.e., as good as other 
models) correlation at the tensor level but give zero solenoidal force, hence zero 
effect on the dynamics of the filtered velocity field, e.g., 

Ti/ 1 - C A 2 (“*“*) • 

Such models should, of course, never be used. Nevertheless, they make 1 the point 
that correlation at the tensor level does not necessarily mean good dynamics. What 
really matters is the solenoidal forcing. 

Finally, scale similarity models of Bardina’s type. Models B, were also investigated 
(e.g., see Horiuti 1993, Zang et al. 1993, Liu et al. 1994, Sal vet ti &: Banerjee 1995). 
The combinations considered were: 


c 

( — 

\ 

(«) 

( UiUj 

-UiUj] , 

c 

[UiUj 

-«|Wj) , 

(b) 

c 

(uiUj 

-«.•«>) . 

(c) 

c 

(ti,Uj 

-UiUj) , 

(d) 


where the additional filtering of the LES field was done at twice the size of the 
original filter, and with the same filter type. Model Ba is the filtering of u t Uj -- T7 t 77 J 
and is ‘similar’ to r t j y which is the original filtering of UiUj — u t TIj. Model Bb is 
similar to r tJ = u t Uj — UjUj. Model Be was also investigated, in a prion testing, 
in Meneveau & Lund (1992). Model Bd is the remainder, after second filtering, of 
u t Uj, and is thus not expected to perform well. 


2. Optimization of the coefficient (s) and correlations 

We first consider optimization of the model’s coefficients ) when no spatial vari- 
ation is allowed. This is a natural requirement as the models are being compared 
with DNS data in isotropic turbulence. For each level, the value of C that pro- 
vides the minimum error, in the least square sense, between ‘exact’ and ‘modeled’ 
subgrid-scale quantities is evaluated. Defining 


_* M 
O 


Cr m-. 


U ’ 


9i 


i cM 


Cg 711 J , 


and (■) as integration over physical space, one then obtains: 
At the traceless tensor level: 


(jh mji) 

( m ti m ki ) 
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At the solenoidal force level: 


C a = 


(nmt) 


' 9 (ml ml) ' 

At the scalar level (= dissipation level): 


d = 




u k9k 


d = u k g k M =Cdm* , 


(d m*) 

d (m*in*) 


For models with two coefficients, defining 


A/ 


C T i + C r2 , 


the least-square optimization then leads to 


( ("»«”*«> ( m h n ki)\ ( C n\ = ( ( T h m ki) \ 

V (™kl n kl) ( n kl n kl) ) \Cr2 / \(rl,n* kl ) ) ' 

and similarly for the coefficients at the force and dissipation levels. 

The correlations between 'exact’ and 'modeled’ subgrid-scale quantities are de- 
fined in the usual way, 


Vr = 


7 h = 


Id = 


( T kt r kt M ) 

m i M ) 

(r k m) ll2 (rk M r k M ) 112 ' 

(T?“) 

(?T) i/2 (d' M d' M )' /2 ’ 


One interesting question is posed when using mixed models instead of the simplest 
models (e.g., using Model 1+3 instead of Model 1 alone or Model 3 alone): by how 
much can one expect to improve the correlation? Part of the answer lies in the 
following identity: 


(q m+n ) 2 = + (JiJl !ll T Jj» 

i-(C ) 2 

where q b stands for the correlation between a and b. Upper and lower bounds are 
obtained as: 


max ((I l ?) 2 M n r?) < (v"r ,+n ) 2 <(v " 1 ) 2 + (’fr ) 2 - 

The above formulas also hold for correlations at the force or energy dissipation levels 
by replacing r/ T by rjg or i ]d- The best situation occurs when the two terms used 
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in the models (a) each correlate well with the exact subgrid-scale quantity, and (b) 
are not highly correlated with each other. This is unfortunately not the case for all 
models tested in this paper (see results), as well as in other papers (e.g., Lund & 
Novikov 1992). 

We then go on to consider local optimization of the model's coefficient s) where 
spatial variation is allowed. Although one usually uses dynamic LES models with 
coefficients ) that, are averaged (and hence uniform) in the directions of flow ho- 
mogeneity (here all three directions), the present study is justified by the hope 
that some model, in some formulation, might exhibit a better behavior than the 
others in terms of coefficients) uniformity. In the same spirit, it is believed that 
dynamic models will then have to work less hard, e.g., require less averaging of the 
coefficients) obtained dynamically. 

The above least square optimization is then carried out locally, at the tensor level. 
For models with one coefficient, the local optimization leads to 

^ _ T h m kt 

^ r — — * — * • 

m kl m kl 

For models with two coefficients, the linear system is solved locally to determine 
Cn and C r2 . The local optimization can be done only at the tensor level. Indeed, 
the force, being the divergence of the sgs tensor, must remain written in conser- 
vative form (so as to have zero global integral in the case of isotropic turbulence). 
Since C r now depends on space, it cannot be pulled out of a derivative such as 
d ( C T m*j ) jdxj. (If one were to relax the constraint that the force be conservative, 
an optimization at the force level would, of course, be possible.) 

Correlations are then computed in the same way as above, at all three levels. At 
the tensor level, they are now artificially much higher than those obtained with uni- 
form C . Since the C s are optimized locally, their spatial variation is quite high, in- 
cluding regions of negative values. Interesting quantities are evaluated and reported: 
mean, rms ((C 2 ) — (C) 2 ) , and normalized pdf, the better models exhibiting a 

sharper (i.e. smaller rms) pdf which is also more skewed towards positive C . 


3. Vorticity-velocity formulation and models investigated 

Consider LES in the vorticity-velocity formulation: 

d 

. -1- I/. J >i — i. 5 ?/ 1 -J- =- =Z 1/ - 

df 


du>, d r 


7 dxj -'dxjdxj 
with the subgrid-scale antisymmetric tensor defined as 


where 


def 



1 (u7|Uj ojjiii) — C t j T Tl t j T L tJ 

7 «> = 

_ def 

7t; “ 

(u >! ,Uj 

-UJjU,) 

( U/ 1 U j 

\jS j U | j 

1 (uJj't/j ~ C i j T j , 

\ / - \ 

C,j = 

("i U 'j 

- “X 

J + ( u l\Tlj - Lj'jUi) , 

R ij = 

*y } - 

-u'ju'i , 

| ( UJ i U j W j ll i ) — ( LiJ i Uj jU t ) 

II 

(uJjjl/j 

u?j U | 
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Scaling 

C T 



warn 




(Cr) 








a 




mm 



0.34 

0.014 





b 





mu 


0.32 

0.016 





c 






0.28 

0.33 

0.23* 


0.57 









0.20 

0.27 

0.012 



Ijlgjl 


HSU 





0.10 

0.26 

0.32 

0.036 

5.8 




d + 








0.032 












0.27 

0.0082 









SIS 



0.31 

0.028 










mm 


0.33 

0.023 

MEM 





Table 1. Investigation of the influence of scaling. Model 1, 128 3 -+ 32 3 with sharp 
cutoff. (*: value of Cf 1 ^ 3 ). 


LES in the vorticity-velocity formulation is a natural choice which requires mod- 
eling only three quantities. Defining f3, = q 23 = -y 32 , /? 2 = 7 3 , = -7^ 
£3 = 712 = -721, together with 7 n = q 22 = 7 33 = 0, one obtains V • 7 = V x 0. 
Hence, modeling is already at the ‘vector’ level, since modeling 7 is really modeling 
Moreover, it is already in the form of a solenoidal ‘forcing’, V x /?. For instance, 
the equivalent of the classical Smagorinsky’s model is simply: 

Model 1 (eddy-viscosity type): 


lij M — —2 Vt r%j with 
which is identical to doing: 


1 fduj, duj\ 

r,j ~ 2 ~ fo~) ’ 


du> 


+ V • (u>u — uuj) = V 








Another model investigated is: 
Model 2 (hyper eddy-viscosity type): 


-A 2 V x (z7, V 4 u) . 


4. Results and discussion 

The velocity-pressure LES formulation is considered first, with the classical Model 1 
but different scalings. The results are compiled in Table 1 for the case 128 3 — ♦ 32 3 
with sharp cutoff. 

For models that can produce local negative scaling, s, it is always better to restrict 
it to positive values by using |s| or s+ = max(s,0): better correlations are obtained 
in the case of uniform C, and sharper distribution of C in the case of variable C T 
(smaller ratio Tms/mean’, e.g., 4-6 instead of 15-16). Of the two options, |s| is 
always slightly better than s + . 

For cases with uniform C (optimized at each level), the correlations are different 
at each level, typically t] T ss 0.12, r] g ss 0.27, and r/j % 0.33. The choice for the 1/T 
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scaling is found to be unimportant as it does not affect the correlations significantly. 
Since the model represents only a crude estimate 1 of the sgs stresses, the correlations 
obtained at the tensor level are quite low. Things improve a bit when considering 
correlation at the force level, and even more so at the dissipation level. This is 
to be expected as one has to work less hard in the modeling effort: 5 degrees of 
freedom versus 3 versus 1. Results in Lund & Novikov (1992) and Clark et al 
(1979) for isotropic turbulence report g T ^ 0.2 (such correlations are also obtained 
here, as presented at the end of this section), Piomelli et al. (198S) for turbulent 
channel flow, and McMillan &: Ferziger ( 1979) for homogeneous shear flow reporting 
rj T ^0.1. In the cast 1 of uniform C\ the correlations only measure the ‘alignment 
between the model and the exact quantities, not the magnitude. Indeed, C drops 
out of the three equations defining the correlations. In particular, one couldmise a 

value which is such that the global dissipation obtained with the model, (d ), be 
equal to the exact global dissipation, (d ): C d — (d These values are also 

reported in the tables. 

For cases with variable C T (optimized locally at the tensor level), the correlations 
are pretty much the same at all levels. To attain correlations in the range 0.56-0.60, 
the model has to ‘work hard 1 : highly varying C r field as seen in the ratio rms/mean 
4-5 and in the normalized pdf of Fig. 3. 

An investigation is also carried out to evaluate the relative participation of the 
two terms C* ? and i?* in in the correlations, see Table 1 2. As expected from the 

mathematical definitions of these terms, the model correlate's better with C l} than 
with R* y It also correlates better with C*j + R l} than with C tJ alone. 


4 

4 

Vr 

Vg 

m 

no 

yes 

0.060 

0.102 

0.174 

yes 

no 

0.108 

0.278 

0.317 

yes 

yes 

0.116 

0.286 

0.339 


Table 2. Contributions of C*j and R*j to the correlations; Model 1 with scaling (a); 
128 3 — > 32 3 with sharp cutoff. 


Model 2 is considered next. This model is not obtained by taking the curl of 
Model 1 (it would if Vt were uniform). Nevertheless, since the choice of scaling was 
found to be unimportant with Model 1, which means that very local variation of v t 
are unimportant, Model 2 is expected to perform as Model 1. This is indeed found 
to be the case. For scaling (a): C g — 0.020, C\i = 0.016, i] g = 0.29, gj = 0.34. 

Model 3 is considered next. The results are reported in Table 3 and in Fig. 3. 
This model performs significantly better than Model 1 at all three levels, and for 
both cases of uniform and variable coefficient. In particular, with variable C r , the 
pdf is sharper (ratio rms/mean of 3.3 instead of 4.4) and more skewed to the right. 

Model 1+3 with uniform coefficients dot's not perform substantially better than 
Model 3 alone (g g — 0.344 instead of 0.334). This is due to the fact that S t } 
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c/c n 


Figure 3. Normalized pdf for Model 1 (• ) and Model 3 (x ) with scaling (a); 
12S 3 — + 32 3 with sharp cutoff. 


Model 

C T 

c g 

c d 

Vr 

Vg 

V d 

(Cr) 

rins 


*19 

Vd 

mean 

Vt 

1 

0.0086 

0.016 

0.013 

0.12 

0.29 

0.34 

0.014 

4.4 

0.57 

0.60 

0.59 

3 

0.0030 

0.0037 

0.0034 

0.16 

0.33 

0.39 

0.0040 

3.3 

0.59 

0.62 

0.63 

14-3 

-0.0077 

-0.014 

-0.012 

0.17 

0.34 

0.40 

-0.0059 

24 

0.75 

0.76 

0.76 


0.0046 

0.0062 j 

i 

0.0060 




0.0053 

5.6 




4 

0.0089 

-0.073 

-0.076 

0.031 

-0.21 

-0.29 

0.019 

12 

0.57 

0.54 

0.53 

1+4 

0.014 

0.015 

0.012 

0.16 

0.29 

0.34 

0.019 

4.6 

0.75 

0.74 

0.73 


0.039 

-0.011 

-0.012 

! 




0.047 

6.9 





Table 3. Comparison of Model 3 and Model 1+3, Model 4 and Model 1+4, with 
Model 1; scaling (a); 128 3 — * 32 3 with sharp cutoff. 


and V 2 5 1; are highly correlated to each other (relative correlations t) T = 0.87, 
T] g = 0.94, rjd = 0.95). It is found that this model corresponds to diffusion with 
the hyperviscosity term (C 2 > 0) corrected by some antidiffusion with the viscosity 
term (C\ < 0). In the case with variable coefficients, it is found that C Tl must vary 
a lot in order to bring the correlations to 0.75 0.7G: very high ratio of rms/mean 
for C r 1 • 

Finally, Model 1+4 is investigated. Model 4 alone performs very poorly: in the 
case of uniform coefficient, the optimization leads to a coefficient of different sign 
whether the correlation is formed at the tensor level, or at the force and dissipation 
levels; in the case with variable coefficient, a ratio rms/mean of 12 is obtained. 
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Model 1+4 with uniform coefficients does not perform better than Model 1 alone. 
This is in accordance with results in Lund k Novikov (1992). This is due not only 
to the high correlations between the two terms (relative correlations 7 ] r = —0.55, 
r) g = —0.67, rjd = —0.82), but also to the low correlations of Model 4 alone. In that 
respect, Model 1+3 performs significantly better than Model 1+4. In the case with 
variable coefficients, the correlations for Model 1+4 axe in the 0.73-0.75 level, with 
not much of an increase in the ratio rms/mean for C rl , and even a decrease in the 
ratio for C r2 • This is in contrast with what was obtained with Model 1+3. 

The case where smooth filters are applied to the DNS data is also investigated. 
In that case, correlations are obtained at the tensor level only. The reason being 
that the filtered data for the r x j was computed on the original 128 3 grid, but was 
only sampled, for output, on a 32 3 subset of that grid. Since it still contains 
very significant contributions from all original modes, this data cannot be properly 
differentiated to obtain the forces. All possible contributions to the sgs tensor were 
considered and were correlated with Model 1, see Table 4. 



c * 

% 

c 

Vr 

(Cr) 

rms 

It 




mean 

no 

no 

yes 

0.00062 

0.083 

0.0014 

7.2 

0.41 

no 

yes 

no 

0.0054 

0.080 

0.011 

12 

0.42 

no 

yes 

yes 

0.0060 

0.084 

0.012 

11 

0.42 

yes 

no 

no 

0.0047 

0.067 

-0.0027 

47 

0.41 

yes 

no 

yes 

0.0053 

0.077 

-0.0014 

92 

0.41 

yes 

yes 

no 

0.010 

0.27 

0.0083 

2.8 

0.34 

yes 

yes 

yes 

0.011 

0.26 

0.0097 

2.6 

0.33 


Table 4. Gaussian filtering of the 128 3 DNS; Model 1 with scaling (a). 


We feel that the cutoff filter is the most appropriate for generating synthetic LES 
fields from DNS data as it completely eliminates the ‘small scale’ information that 
will never be present in a large eddy simulation. In that sense, a priori tests using 
smooth filters such as the Gaussian are of a more academic interest. Indeed, an 
LES simulation would not be able to capture with the 32 3 grid the ‘small scale’ 
information which is still present after smooth filtering of the 128 3 DNS data. With 
the Gaussian, the filter value at the cutoff wavenumber (i.e., at the edge of the 32 3 
grid) is 0.663 3 = 0.291, which is still very significant. At twice that wavenumber, it 
has dropped to 0.193 3 = 0.0072. It thus can be argued that a 64 3 grid (or so) would 
be needed to correctly capture the important part of the fine grain information left 
after Gaussian filtering. 

Proceeding nevertheless^ with this study, one finds that the case where Model 1 
is correlated with C* tj + R* tj + L*j (= r£) performs quite well: r) T = 0.26 for the 
case with uniform coefficient, rms/mean of only 2.6 to reach r) T = 0.33 for the 
case with variable coefficient. Similarly for the case C* tJ + L*j f . The cases L*j 

and Rij + L perform very poorly: very low correlation in the case with uniform 
coefficient, extremely high value of rms/mean in the case with variable coefficient. 
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C/C mean 

FIGURE 4. Normalized pdf for Model 1 with scaling (a); Gaussian (• ) and top-hat 
(■ ) filtering of the 128 3 DNS. 

All other case are such that L*j is not included and perform very poorly, e.g., the 

case C i} + R tJ ( = r\ } ) , which has only r; r = 0.084. In conclusion, the correlation 
with r*j is significantly higher than what was obtained with sharp cutoff (t} r = 0.26 
instead of 0.12) but the correlation with is significantly lower (r} T = 0.084 instead 
of 0.12). Finally, very similar results are obtained when using the top-hat filter 
instead of the Gaussian, see e.g., Fig. 4 for very similar pdf’s. 

The scale similarity models, Models B, are considered next, first with sharp cutoff 
filter, then with Gaussian filter. 


Type 

C T 

^9 

c d 

c d ' 

V r 

ng 

Vd 

a 

0.028 

0.022 

0.084 

0.68 

0.043 

0.017 

0.070 

b 

0.015 

0.031 

0.083 

0.43 

0.039 

0.034 

0.098 

c 

-0.0054 

-0.044 

-0.053 

1.20 

-0.022 

-0.12 

-0.14 

d 

0.0076 

0.039 

0.090 

1.2 

0.015 

0.031 

0.075 


Table 5. Scale similarity models, 128 3 — ► 32 3 with sharp cutoff. 

With sharp cutoff, Models B are only compared to t* } (since cannot be 
captured on the 32 3 grid), see Table 5. It is seen that very low levels of correlations 
are obtained regardless of the Model’s type. This is in accordance with results by 
Meneveau &: Lund (1992). Again, one must recall that, in the case of uniform C\ 
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the correlation is independent of C. Nevertheless, when C is optimized as usual, one 
obtains a very large difference between C,/ and C T , C g . Cj. which is also indicative 
of very poor models. Thus, with sharp cutoff filter, Bardina’s type models show 
essentially no correlation with the relevant sgs quantities: trace-free tensor t*j , 
solenoidal force g* . In fact, if anything, Models Ba and Bb only correlate a little 
when the trace of t,j is kept (g T = 0.13 and 0.11 respectively). 


Type 

sgs 

C T 


(Cr) 

rms 


'Hr 

mean 

7 It 

a 

T ij 

0.14 

0.29 

0.21 

1.8 

0.61 


n, 

0.57 

0.70 

0.51 

0.93 

0.73 

b 

Tti 

0.47 

0.70 

0.38 

0.90 

0.74 


n, 

0.29 

0.25 

0.33 

2.8 

0.53 

c 


0.014 

0.062 

0.0073 

48 

0.47 



0.37 

0.92 

0.36 

0.57 

0.78 

d 


0.095 

0.20 

0.15 

2.8 

0.57 



-0.47 

-0.55 

-0.40 

1.44 

0.64 


Table 6 . Scale similarity models; Gaussian filtering applied to the 128 3 DNS. 

With Gaussian smoothing, things are completely different, see Table 6 . Models B 
are here compared to both rf and t*j. Model Ba is ‘similar’ to t*j and correlates 
indeed very well with it: tj t = 0.70 for the case with uniform C, very low rms/mean 
of 0.93 for the case with variable C. Model Bb is ‘similar’ to r* and correlates 
indeed very well with it: ij T = 0.70 for the case with uniform C , rms/mean of 
0.90 for the case with variable C. Model Be correlates very poorly with r,* , but 
extremely well with T* } : g T - - 0.92 for the case with uniform C, rms/mean of 0.57 for 
the case with variable C. This is the highest correlation encountered in the course 
of this study. It is consistent with the 0.8 correlation reported in Meneveau k Lund 
(1992). Model Bd does not perform well, as expected, since it is the remainder of 
some quantity after second filtering. 

These impressive results are misleading. Indeed, Gaussian filtering of the DNS 
data produces a synthetic LES field that still contains considerable contributions 
from the small scales. As this small scale information will not be present in an real 
LES, results obtained with the sharp cutoff filter are more representative of what 
might be expected from using Bardina’s models in LES. 

The vorticity- velocity LES formulation is now considered. The correlation at the 
antisymmetric tensor level, 7 is denoted as C 7 , at the ‘forcing level as C g . at the 
‘enstrophy dissipation’ level as C',{ . The results obtained with the eddy- viscosity 
model, Model 1, and with the hyper eddy- viscosity model, Model 2, are presented 
in Table 7. 

The eddy-viscosity model in the vorticity-velocity formulation produces signifi- 
cantly higher correlations than its counterpart in the velocity-pressure formulation: 
in the case of uniform coefficient, rj-, = 0.19 instead of r/ T = 0.12, r) g = 0.32 instead 
of 0.29, ijd = 0.46 instead of 0.34; in the case of variable coefficient, = 0.71 
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Model 

c. 

c. 

c , 

CJ 

Vy 

Vg 


(Cy) 



Vg 

Vd 

U d 


1 

0.019 

0.026 

0.019 

0.021 

0.19 

0.32 



4.6 

0.71 

0.74 

0.80 

2 

0.0041 

0.0046 

0.036 

0.039 

0.23 

0.35 

0.48 

0.0045 

4.0 

0.73 

0.75 

0.79 


Table 7. Vorticity-velocity formulation; Model 1 and Model 2 with scaling (a); 
12S 8 — > 32 3 with sharp cutoff. 


instead of r; r = 0.57, i ] g = 0.74 instead of 0.60, r)d = 0.80 instead of 0.59, with 
essentially the same ratio rms/mean as before, and a pdf which is more skewed 
towards positive C\ see Fig. 5. 

Again, the hyper eddy-viscosity version of the model performs even better than 
the eddy-viscositv version, see Table 7 and pdf of Fig. 5. 



FIGURE 5. Normalized pdf for Model 1 (■ ) and Model 2 (+ ) with scaling (a) in 
the vorticity-velocity formulation; 128 3 — > 32 3 with sharp cutoff. 

Moreover, one finds for both models that the coefficients optimized globally are 
close to each other, and that they are also close to the average of the coefficient 
optimized locally. This is also indicative of good candidate models for LES. 

For completeness, the case of smooth filtering of the DNS data is also investigated, 
see Table 8. The Gaussian and top-hat filter produce similar results. The eddy- 
viscosity model in the vorticity-velocity formulation performs slightly better than 
in the velocity-pressure formulation: in the case of uniform coefficient, r\ 1 = 0.28 
instead of 0.26; in the case of variable coefficient, r/ 7 = 0.55 instead of 0.33, with 
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Type 

sgs 

C\ 

r h 

(C-r) 

rms 

>h 

mean 

Gaussian 

li] 

0.014 

0.28 

0.014 

2.1 

0.55 

top-hat 

Hj 

0.015 

0.24 

0.016 

2.3 

0.53 

Gaussian 


0.042 

0.051 

0.0009 

140 

0.52 

tophat 

_lii_ 

0.050 

0.057 

0.0029 

42 

0.53 


Table 8 . Vorticity- velocity formulation; Gaussian and top-hat filtering of the 128 3 
DNS; Model 1 with sealing (a). 


a smaller ratio rms/mean (2.1 instead of 2.6). Again, when the L' tJ term is not 
included (i.e., when using instead of 7 ^), the model correlates very poorly. 

So far, all correlations have been obtained using the same 128 3 DNS database in 
forced isotropic turbulence (Re\ = 95.8). With sharp cutoff filter, the filtering was 
always done from 128 3 DNS to a 32 3 synthetic LES field. An investigation is now 
done (a) using the same database, but filtering to a 16 3 synthetic LES field, and (b) 
filtering the 512 3 DNS database in decaying isotropic turbulence (Re\ = 63.5) to 
both 64 3 and 32 3 synthetic LES fields. The regions of the spectra where the cut is 
done are marked in Figs. 1 and 2 . Although there is no pure A * -5 / 3 'inertial range’ in 
this 512 3 computation at such Re\, there is 'almost’ an inertial range, the 64 3 cut 
being to the far right of it, and the 32 3 cut being within it. With the 128 3 database, 
the 32 3 cut is also to the far right of the 'inertial range’ (actually probably more 
at the beginning of the 'dissipation range’), and the 16 3 cut is within the inertial 
range. Results of this investigation done in the velocity-pressure formulation are 
reported in Table 9 and in Fig. 6 . The results corresponding to the vorticity- velocity 
formulation are presented in Table 10. 


Data 

cv 

Cg 

c d 

Cj 

Vr 

Vg 

m 

(CV) 

rms 


Vg 

rid 

mean 

r h 

128 3 -► 32 3 

0.0086 

0.016 

0.013 

0.0098 

0.12 

0.29 

0.34 

0.014 

4.4 

0.57 

0.60 

0.59 

128 3 -> 16 3 

0.061 

0.077 

0.070 

0.069 

0.18 

0.32 

0.46 

0.11 

3.5 

0.49 

0.38 

0.38 

512 3 -» 64 3 

0.016 

0.019 

0.018 

0.018 

0.18 

0.30 

0.42 

0.022 

3.1 

0.59 

0.60 

0.64 

512 s — > 32 * 

0.029 

0.028 

0.028 

0.031 

0.29 

0.39 

0.59 

0.036 

2.0 

0.61 

0.64 

0.73 


Table 9. Investigation of different databases and of different cut locations in each 
database; velocity- pressure formulation; Model 1 with scaling (a); sharp cutoff. 


We concentrate on correlations obtained with uniform coefficients. When con- 
sidering different cuts within the same database, one finds that the cut within the 
inertial range produces higher correlations than the cut to the far right of that 
range. We don’t see any obvious reason at this time why this should be the case. 
Nevertheless, this finding holds for both databases investigated and for both formu- 
lations. 

Notice that the superior performance of the vorticity- velocity formulation over 
the velocity-pressure formulation is not as marked in the 512 3 runs as it is in the 
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Data 


c 


c V 


Vg 

f ld 

(C,) 

rms 

Vt 

J lg 

Vd 

^ 9 

mean 

128* -+ 32 3 

0.019 

0.026 

0.019 

0.021 

0.19 

0.32 

0.46 

0.023 

4.6 

0.71 

0.74 

0.80 

512 3 — 64 3 

0.022 

0.025 

0.023 

0.024 

0.21 

0.29 

0.50 

0.025 

4.3 

0.71 

0.73 

0.81 

512 3 — 32 3 

0.032 

0.033 

0.033 

0.035 

0.28 

0.35 

0.62 

0.037 

3.0 

0.73 

0.74 

0.84 


Table 10. Investigation of different databases and of different cut locations in each 
database; vorticity-velocity formulation; Model 1 with scaling (a); sharp cutoff. 



^ ^ mean 

FIGURE 6. Normalized pdf for Model 1 with scaling (a) in the velocity-pressure 
formulation; 128 3 -*• 32 3 (v ), 128 3 -> 16 3 (* ), 512 3 -► 64 3 (■ ), 512 3 4 32 3 (• ) 
with sharp cutoff. 

128 3 runs. In particular, the 512 3 — ► 32 3 cases are very comparable when C is 
uniform (still slightly better when C is allowed to vary). 

5. Conclusions 

A few conclusions can be made from the investigations done using the sharp cutoff 
filter to produce the synthetic LES fields from the DNS databases. The choice for 
the 1/T scaling in the eddy- viscosity is found to be unimportant as it does not 
significantly affect the? correlations between modeled and exact sgs quantities at any 
of the three levels: trace-free tensor, solenoidal forcing, dissipation. It is found that 
the hyper eddy-viscosity model yields higher correlations than the eddy-viscosity 
model in both the velocity-pressure and the vorticity-velocity formulations. It thus 
appears as a good candidate for real LES and should be tested numerically. Scale 
similarity models exhibit essentially no correlation with the exact sgs quantities 
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when the sharp cutoff filter is used. The correlations obtained with simple LES 
models in the vorticity-velocity formulation are higher than those obtained with 
the counterpart models in the velocity-pressure formulation. This suggests that the 
vorticity-velocity formulation might be a good candidate for real LES, with reduced 
sgs modeling error. It certainly appears as a natural choice. It should also be tested 
numerically in real LES. 

Some conclusions are also reached from the investigations done using the Gaus- 
sian and top-hat filters in order to produce the synthetic LES fields. It is found 
that the L' t] contribution is essential in order to produce significant correlations. 
The correlations are then artificially higher than those obtained with sharp cutoff. 
Indeed, the filtered field still contains significant small-scale information that would 
not be available in a real LES. In particular, some scale similarity models exhibit 
remarkably high correlations when smooth filters are used. 

Finally, it is found that the level of obtained correlation is quite sensitive to the 
database investigated, and to the location of the spectral ‘cut’ used to produce the 
synthetic LES field. If anything, our investigation shows that one must exercise cau- 
tion when comparing correlations reported by different authors and when working 
with different databases and with different models or formulations. 
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